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ABSTRACT 

Recent observations of high-redshift galaxies and quasars indicate that the hydrogen 
component of the intergalactic medium (IGM) must have been reionized at some redshift 
z ^ 6. Prior to complete reionization, sources of ultraviolet radiation will be seen behind 
intervening gas that is still neutral, and their spectra should show the red damping wing 
of the Gunn-Peterson trough. While this characteristic feature may, in principle, totally 
suppress the Lya emission line in the spectra of the first generation of objects in the 
universe, we show here that the IGM in the vicinity of luminous quasars will be highly 
photoionized on several Mpc scales due to the source emission of Lyman-continuum 
photons. If the quasar lifetime is shorter than the expansion and gas recombination 
timescales, the volume ionized will be proportional to the total number of photons 
produced above 13.6 eV: the effect of this local photoionization is to greatly reduce the 
scattering opacity between the redshift of the quasar and the boundary of its H n region. 
We find that the transmission on the red side of the Lya resonance is always greater 
than 50% for sources radiating a total of > io 69 - 5 ionizing photons into the IGM. The 
detection of a strong Lya emission line in the spectra of bright QSOs shining for > 10 7 yr 
cannot then be used, by itself, as a constraint on the reionization epoch. The first signs 
of an object radiating prior to the transition from a neutral to an ionized universe may 
be best searched for in the spectra of luminous sources with a small escape fraction of 
Lyman-continuum photons into the IGM, or sources with a short duty cycle. 

Subject headings: cosmology: theory - galaxies: formation - intergalactic medium - 
quasars: absorption lines - radiative transfer 

1. Introduction 

Popular cosmological models predict that most of the intergalactic hydrogen was reionized by 
the first generation of stars or accreting black holes in the universe at a redshift 7 < z re \ < 15 
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(e.g. Ciardi et al. 2000; Haiman & Loeb 1997; Gnedin & Ostriker 1997). Studying the epoch of 
reheating which marks the end of the cosmic "dark ages" is crucial for determining its impact on 
several key cosmological issues, from the role reionization plays in allowing pregalactic objects to 
cool and make stars, to determining the small-scale structure in the temperature fluctuations of the 
cosmic microwave background. Because of scattering off the line-of-sight due to the diffuse neutral 
intergalactic medium (IGM), the spectrum of a source at z em > z re \ should show a Gunn-Peterson 
(1965, hereafter GP) absorption trough at wavelengths shorter than the local Lya resonance, A Q bs < 
A«(l + 4m)i where \ a — cjv a — 1216 redA. Unfortunately, while actively being searched for, 
the first signs of an object radiating prior to the reionization epoch are far from unambiguous. 
This is because: (1) line blanketing from discrete Lya forest absorbers becomes really severe at 
z > 5.5, with less than 10% of the quasar unabsorbed continuum leaking through between rest- 
frame Lya and Ly/3 (Fan et al. 2000a; Stern et al. 2000); and (2) the GP optical depth for resonant 
scattering encountered by photons propagating through smoothly distributed neutral hydrogen of 
density nm(z), 

blue/ n ™hi /x vre 2 /A a _ inSfc-irv-Va f ^bh 2 \ /l + z\ 3/2 

r GP (z) = r — , ^) = __ nH= , ujxiofc n M ' {—) , (1) 

is extremely high. Here / is the oscillator strength, e and m e are the electron charge and mass, % 
is the mean density of hydrogen nuclei at redshift z, Hq = 100 h kms -1 Mpc -1 is the present-day 
Hubble constant, and (Oj^jfij) are the total matter and baryonic density parameter, respectively. 
The expression above gives the opacity seen by any photon emitted from a source at z cm on 
the blue side of the Lya line, A em < A a , as it is redshifted through the local Lya resonance 
at (1 + z) = (1 + -Zcm)A C m/A Q . Equation (|l|) shows that the transmitted quasar flux shortward 
of Lya would be reduced to undetectable levels even if 99% of all the cosmic baryons were to 
fragment at these early epochs into discrete, mildly overdense structures, with only 1% remaining 
in a diffuse component that was 99% ionized. As such, the detection of a GP trough would not 
uniquely establish that an object is being observed prior to the reionization epoch (except, perhaps, 
in the case where reionization occurs extremely rapidly and the GP trough splits into individual 
Lyman series troughs for a source located at (1+^rei) < (l + ^em) < 32(l + 2: re i)/27, Haiman & Loeb 
1999). It has been pointed out by Miralda-Escude (1998), however, that the rest-frame ultraviolet 
spectra of sources observed prior to complete reionization should show the red damping wing of 
the GP trough, as they will be seen behind a large column density of intervening gas that is still 
neutral. At z > 6, this characteristic feature extends for more than 1500 kms -1 to the red of the 
resonance, significantly suppressing the Lya emission line. Measuring the shape of the absorption 
profile of the damping wing could provide a determination of the density of the neutral IGM near 
the source. 

In this Letter we focus on the width of the red damping wing - related to the expected strength 
of the Lya emission line - in the spectra of very distant quasars as a flag of the observation of 
the IGM before reionization. We discuss, in particular, the impact of the photoionized, Mpc-size 
regions which will surround individual luminous sources of UV radiation on the transmission of 
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photons redward of the Lyct resonance, and show that the damping wing of the GP trough may 
nearly completely disappear because of the lack of neutral hydrogen in the vicinity of a bright 
object. The absence in the spectra of luminous quasars of a red damping wing with the predicted 
absorption profile will not provide then unambiguous evidence of the observation of the IGM after 
cosmological reionization. 



2. The red damping wing 



We generalize here the calculation of the absorption profile of the damping wing of the GP 
trough (Miralda-Escude 1998) to the case where the IGM in the vicinity of an object at z em is pho- 
toionized due to the source emission of UV photons. When the column density of absorbing atoms 
is sufficiently large, the width of an absorption line substantially exceeds the value corresponding to 
the dispersion of particle velocities along the line of sight. In this case, the scattering cross-section 
is determined by the natural width of the Lya resonance, 
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(Peebles 1993), where A = (87r 2 e 2 /)/(3m e cA 2 f ) = 6.25 x 10 8 s _1 is the rate of spontaneous radiative 
decay from the 2p to Is energy level. The regions of the line profile for which equation (Q) is valid 
are known as "radiation damping" wings of the line. Equation (|2]) is unapplicable close to line 
center, but if damping wings are present the transmitted flux will be essentially zero in the core 
region anyway. We assume the IGM has a constant H I comoving density nni(O) at all redshifts 
z re i < z < Zi < z ern , and is highly photoionized between the redshift of the source z em and the 
boundary of its H II region at Zj. 

The scattering optical depth at the observed wavelength A Q b s > A a (l + z cm ) is 
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where dljdz = c[(l + z)H(z)] 1 is the proper cosmological line element. In an EdS universe, 
equation (||) can be rewritten as 
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where R = AAq,/(47tc) = 2.02 x 10~ 8 , x re i = (1 + 2 re i)AQ,/A b s , and X{ = (1 + Zi)\ a / \ -^ s . Far from 
line center (i.e. when || 1 — x ||^> Rx 3 ), this integral has an analytic solution (Miralda-Escude 
1998). Figure 1 (solid curve) shows the red damping wing in the spectrum of a source at z em = 7 
assuming T (z cm ) = 3 x 10 , 

Zrei — 6, and zi — z em , i.e. in the case the H n region surrounding 
the radiation object is very small (because, e.g., the Lyman-continuum photons produced cannot 
escape from the dense sites of star formation into the intergalactic space). 
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Equation (||) also gives the opacity at wavelengths X a (l + z\) < A D bs < A Q (1 + z cm ), i-e. on the 
blue side of the quasar Lya emission line, due to the damping wing of the fully neutral gas along 
the line of sight at z roi < z < z\. We will see in the next section that this must be augmented by 
the scattering optical depth of the residual H I in the vicinity of the source (zj < z < z em ). 



3. Cosmological H II regions around isolated sources 

We now assess in details the impact of a local H n region on the shape of the damping wing 
profile. When an isolated point source of ionizing radiation turns on, the volume of ionized IGM 
initially grows in size at a rate fixed by the emission of UV photons, and an ionization front 
separating the H n and H I regions propagates into the neutral gas. Most photons travel freely in 
the ionized bubble, and are absorbed in a transition layer (the "I-front"), across which the degree 
of ionization changes sharply on a distance which is small compared to the radius of the ionized 
zone (this is true even in the case of a QSO with a hard spectrum, see Madau & Meiksin 1991). 
The evolution of an expanding cosmological H n region is governed by the equation 

W-iHr,-*-*, (5) 
at nn t Tec 

(Shapiro & Giroux 1987; Madau, Haardt, & Rees 1999), where Vi is the proper photoionized 
volume, Ni is the number of H-ionizing photons emitted by the central source per unit time that 
escape into the IGM, 

-l 

x 0.019 

is the volume-averaged recombination time, as is the radiative recombination coefficient to the 
excited states of hydrogen (at an assumed gas temperature of 10 K), and the factor C = (rip) /rip > 
1 takes into account the degree of dumpiness of the photoionized region, with n p the mean proton 
density. When the source lifetime t s is much less than (t Tec , H~ l ), as expected for example in 
the case of a quasar shining for an Eddington time scale, t s = tE = 4 x 10 (e/0.1)yr with e the 
radiative accretion efficiency, recombinations can be neglected (this would be true even in the case 
of a clumpy medium with C < 10 on Mpc scales, cf eqs. || and ^) and the evolution of the H n region 
can be decoupled from the expansion of the universe. The volume, 

(t s «t rcc ,^ 1 ), (7) 

that is actually ionized becomes then proportional to the total number of Lyman-continuum photons 
emitted over the source lifetime. Note that the I-front initially expands at velocities that are close 
to the speed of light, i.e. with rj(t) = [3V/(t)/47r] 1//3 ~ ct. In the case of very luminous, short-lived 
objects this 'relativistic' phase may last a considerable fraction of the source lifetime. Nevertheless, 
the apparent size of the H n region 'seen' by Lya photons propagating along the line of sight will 
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always be given by equation (^), as these can catch up with the I- front only when it slows down to 
subluminal velocities. 

The known QSOs at z em > 5 all have ionizing luminosities that can be estimated to lie in the 
range 10 56 < iVj < 10 58 s^ 1 , the brightest of them being the recently discovered z em = 5.8 quasar 
from the Sloan Digital Sky Survey (SDSS, Fan et al. 2000a). Over a lifetime of (say) 10 7 yr they 
would radiate of order N{ = 10 70 ' 5 — 10 72 ' 5 photons above 13.6 eV. Placed at z em = 7, sources 
of similar power would ionize the surrounding IGM out to a proper distance 77 = 1.5 — 7Mpc, 
corresponding to a Hubble expansion velocity of Av = Hrj = 3300 — 15800 hQ, 1 ^ 2 kms -1 , or to a 
redshift difference between the QSO and the boundary of its H n zone of Az = 3.33 x 10 -4 h r2^ 2 (l+ 
£em) 5/2 (r//Mpc) = 0.088 - 0.42/if^ 2 . The effect of these individual H n regions on the flux 
transmission redward of the Lya line is shown in Figure 1. The width of the damping wing measures 
the column density along the line of sight of the fraction of the IGM that is still neutral. The wing 
nearly completely disappears in the case of a luminous quasar, as there is very little neutral gas in 
its vicinity. In particular, the transmission is always greater than 50% for iVj > 10 69 - 5 photons. 

On the blue side of Lya , in the vicinity of the source, the scattering opacity has two contri- 
butions, one due to the damping wing of the fully neutral gas beyond the H n region, and one due 
to the ordinary GP trough associated with the residual H I in the photoionized zone. Assuming 
photoionization equilibrium (this is justified for the highly ionized IGM and the source lifetimes 
considered here), at every point within the H n bubble the density of neutral hydrogen is given by 
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where hp is the Planck's constant, F u is the incident ionizing flux per unit frequency - to a first 
approximation simply the radiation emitted by the quasar reduced by geometrical dilution - and 
the hydrogen photoionization cross-section (by photons above the threshold hpfi = 13.6 eV) is 

aa{y) ~ a L (v/v L )~ 3 , a L = 6.3 x 10~ 18 cm 2 . (9) 

For a power-law spectrum of the form F u oc v~ a near the hydrogen Lyman edge, the neutral 
hydrogen column through an (approximately isothermal) H II region can thus be written as 

iVm= r drn m = ^a- L 1 ^, (10) 
Jo ol t rec 

independently of the source luminosity. Taking a = 0.5, t s = 10 7 yr, and i rec = 1.3 Gyr at z cm = 7, 
one derives Nm = 8.6 x 10 15 cm~ 2 . This is too small a column (even in the case of a fast recombining 
clumpy medium) for the natural width of the line to exceed the thermal broadening. On the other 
hand, it is straightforward to derive from equations ([[]), ©, and @ that the GP optical depth on 
the blue side of Lya due to partially ionized gas at a distance r < rj from the QSO is given by 
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While the first two terms are of course smaller than unity, the third term is rather large, cH 1 /rj 

r 1 



130 h 1 fl A j^ 2 /rj Mpc at z em = 7, and this means that only the inner parts of the H n region will be 



optically thin to the classical GP absorption. Note that r^p 6 just depends on the quasar luminosity. 



not on t s and rj separately, and that equations (g), (PI, and (11) assume: (1) a steady luminosity, 
and (2) a uniform medium at the mean background density. The hydrogen neutral fraction will 
actually depend on the mean ionizing flux over the last t Tec (ri}ii/np) <C t vec vr, while halos having 
collapsed from (say) 3a fluctuations may actually be sitting in slightly overdense regions. 



4. Summary 

The lack of a GP trough - i.e. the detection of transmitted flux shortward of the Lya wavelength 
- observed in the spectrum of the z em = 5.8 SDSS quasar (Fan et al. 2000a) indicates that the IGM 
was already highly ionized at that redshift, and that sources of ultraviolet photons were present in 
significant numbers when the universe was less than 6% of its current age. At the time of writing, 
three quasars have already been found in the range 5 < z em < 5.5 (Zheng et al. 2000; Stern et al. 
2000; Fan et al. 2000b), and two galaxies have been spectroscopically confirmed at z em > 5.6 (Hu, 
McMahon, &; Cowie 1999; Weymann et al. 1998). It is estimated that the SDSS could reveal tens 
of additional sources at these epochs, and one QSO at z em > 6 about twice as luminous as 3C273 
in every 1500 deg 2 of the survey (Fan et al. 2000a). 

Near-future studies of the rest-frame UV spectra of high-redshift quasars and star-forming 
galaxies could then conceivably provide important probes of the formative early stages of cosmic 
evolution. In this Letter we have discussed one of the first signs of an object radiating prior to 
the transition from a neutral to an ionized universe, the red damping wing of the GP trough. 
We have shown that the local photoionized zones which will inevitably surround luminous quasars 
at early epochs will greatly reduce the scattering opacity between the redshift of the source and 
the boundary of its Mpc-size H n region, thus increasing the transmission of photons on the red 
side of the Lya resonance. To better gauge this effect on real data, we have plotted in Figure 2 
the Keck/LRIS spectrum of the faint z em = 5.5 quasar RD J0301117+002025 (Stern et al. 2000), 
redshifted to z em = 7. The figure depicts the same 800 redA-wide region of the observed spectrum 
around the Lya resonance, together with the transmission exp(— tqp) assuming the QSO is being 
observed prior to the reionization epoch at z re \ = 6. This should be taken just as an illustrative 
example, as in some numerical simulations (e.g. Ciardi et al. 2000) reionization was already well in 
progress prior to redshift 6, and the form of the damping profile would be different in the case of 
patchy ionization along the line of sight. Four cases are shown, as the emission rate of UV photons 
which ionize the IGM in the vicinity of the quasar is increased from TV; = to Ni = 10 56 ,10 57 , 
and 10 58 s . The calculations assume a source lifetime of t s = 10 7 yr. The suppression effect of 
the Lya emission line by the wing of the GP trough, while clearly visible in the absence of a local 
H II region (top-left panel), weakens significantly as the size of the photoionized zone increases. On 
the blue side of the resonance, the transmission profile is due to the combination of two effects: 
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(1) the damping wing - i.e. the curve in the top-left panel at A Q bs > A a (l + z em ) - shifted to the 
blue by the difference between the redshift of the H n region boundary and the quasar redshift; 
and (2) the ordinary GP optical depth from the residual neutral gas in the H n zone. The first of 
these contributions depends on iVj (i.e. on the luminosity iVj multiplied by the source lifetime t s ), 
whereas the second just depends on Ni. Figure 3 shows a similar set of spectra and transmissions 
at earlier epochs, with z em = 9, z vc \ = 8, and a recombination timescale which is (10/8) 3 times 
shorter. 

Depending on the amount of H-ionizing photons which escape the H I galactic layers into the 
intergalactic space, the red damping wing may be more easily detected in the spectra of star- 
forming early galaxies (as well as luminous sources with a short duty cycle). For standard initial 
mass functions, about 3000-4000 photons are produced above 13.6 eV for every stellar baryon. Thus 
a galaxy which turns a gas mass M* into stars will radiate a total of about 10 68,6 (M*/10 8 M Q ) 
Lyman-continuum photons. It is clear from Figure 1 that the effect of its local H II region on the 
shape of the GP absorption profile would be negligible if the escape fraction was (say) < 10%. In 
this case, large H n regions may still be expected if these early galaxies were highly clustered on 
Mpc scales. 

Support for this work was provided by NASA through ATP grant NAG5-4236 (P. M.), by a 
B. Rossi Visiting Fellowship at the Observatory of Arcetri (P. M.), and by the Royal Society (M. 
J. R.). We are indebted to D. Stern and H. Spinrad for providing the spectrum of quasar RD 
J030117+002025. Results similar to the ones presented here have been reached independently by 
Cen & Haiman (2000). 
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Fig. 1. — The red damping wing of the Gunn-Peterson trough for a source at z em = 7. The 
transmission exp(— rg£) is plotted as a function of the fractional wavelength interval from the 
Lya resonance at A Q (1 + z em ). The neutral IGM has tq = 3 x 10 5 at z em (cf. eq. 1), and is 
assumed to be completely reionized by z re \ = 6. Solid curve: absorption profile neglecting the effect 
of the local H n region around the source (Ni = 0). Note how a fraction > exp(— 0.5) = 0.6 of 
the radiated flux will only be transmitted > 1200 kms -1 (corresponding to > 40redA) to the red 
of the resonance. Dashed curves: absorption profiles including the local H n zone generated by a 
QSO shining for 10 7 yr. From top to bottom, the curves are plotted for four decreasing ionizing 
luminosities, iVj = 10 58 , 10 57 , 10 56 , and 10 55 s _1 (Einstein-de Sitter universe with h = 0.5). 
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Fig. 2.— The Keck/LRIS spectrum of the faint z em = 5.5 quasar RD J0301117+002025 (Stern 
et al. 2000), redshifted to 

z em — 7. The dotted vertical lines indicate the expected wavelength of the 
Lya resonance. Fluxes have arbitrary normalizations. The dashed curves show the transmission 
exp(— tqp) through a uniform IGM, assuming the QSO is being observed prior to the reionization 
epoch at z re \ = 6. Four different cases are shown, as the rate of emission of Lyman-continuum 
photons which ionize the IGM in the vicinity of the source is varied from JV f = to N { = 10 56 , 10 57 , 
and 10 58 s _1 . The calculations assume a quasar lifetime of t s = 10 7 yr, a power-law spectrum of 
the form F v oc v~ a with a = 0.5 near the hydrogen Lyman edge, a recombination timescale of 
t rec = 1.3 Gyr, and an Einstein-de Sitter universe with h = 0.5. 
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Fig. 3. — Same as Figure 2, but with z em = 9, z re \ = 8, and a recombination timescale (10/8) 3 
times shorter. 



